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Abstract 

The quasigeostrophic model is a simplified geophysical fluid model 
at asymptotically high rotation rate or at small Rossby number. We 
consider the quasigeostrophic equation with dissipation under ran- 
dom forcing in bounded domains. We show that global unique solu- 
tions exist for appropriate initial data. Unlike the deterministic quasi- 
geostrophic equation whose well-posedness is well-known, there seems 
no rigorous result on global existence and uniqueness of the randomly 
forced quasigeostrophic equation. Our work provides such a rigorous 
result on global existence and uniqueness, under very mild conditions. 
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1 Introduction 



The models for geophysical flows are usually too complicated for analysis. 
Simplified partial differential equation models which are intended to capture 
the key features of large scale phenomena and filter out undesired high fre- 
quency oscillations in geophysical flows have been derived at asymptotically 
high rotation rate or small Rossby number. An important example of such 
partial differential equations is the quasigeostrophic model. 
The deterministic quasigeostrophic equation is ( 150] , B ) 



Aip t + J(ip, Aip) + f3tp x — vA 2 ip — rAtp , (1) 

where ip(x,y,t) is the stream function, (3 > is the meridional gradient of 
the Coriolis parameter, v > is the viscous dissipation constant, r > is 
the Ekman dissipation constant. Moreover, J(f,g) = f x g y — f y g x denotes 
the Jacobian operator. 

The deterministic quasigeostrophic equation (Jj) has been derived as an 
approximation of the rotating shallow water equations by the conventional 
asymptotic expansion in small Rossby number (P0|). Schochet ( P^l ) has 



recently shown that the shallow water flows converge to the quasigeostrophic 
flows in Sobolev norms in the limit of zero Rossby number (i.e., at asymptoti- 
cally high rotation rate), for appropriate initial data. For related work about 
the three dimensional baroclinic quasigeostrophic model, see, for example, 
§, 0, (T7| and §. 

Recently, a few authors have considered the randomly forced quasi- 
geostrophic equation, in order to incorporate the impact of uncertain geo- 
physical forces []13| , fl5fl , fl8|l , fllCfl ). They studied statistical issues 



such as estimating correlation coefficients for the linearized quasigeostrophic 
equation with random forcing. There is also recent work about the impact 
of random ocean bottom topography on quasigeostrophic dynamics (||16||). 
The randomly forced quasigeostrophic equation takes the form 



Aifj t + J(V>, Aip) + Pif> x = uA 2 iJj - rA^j + ^ , (2) 

where W(x, y, t) is a space-time Wiener process to be defined below. There 
does not seem to exist a mathematically rigorous theory of quasigeostrophic 
dynamics under random forcing. In this paper, we consider existence and 
uniqueness of solutions for the nonlinear quasigeostrophic equation (fj) sub- 
ject to Dirichlet boundary conditions and appropriate initial data. 
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2 Local existence and uniqueness of solution 



processes 

Introducing ui(x,y,i) = Aip(x,y,t), the equation (Q) can be written as 

uj t + J(ip, to) + j3ij} x = vAu - ruj + , (3) 

where (x, y) G D and D C M 2 denotes a bounded domain with sufficiently 
regular boundary. This equation is supplemented by zero Dirichlet boundary 
conditions (0) for both if) and u = Aip, together with an appropriate initial 
condition, i.e., we require 

ip(x,y,t) = ondD, (4) 
u(x,y,t) = on dD , (5) 
u(x,y,0) = u (x,y) . (6) 

We note that the Poincare inequality holds with these boundary conditions. 

As it stands, (§) still has to be given a mathematically precise mean- 
ing. This can be done using the framework of stochastic partial differential 
equations (0). For this we (formally) rewrite @ in the form 

du = (vAu -ru- (3ij x - J (if), u))dt + dW . (7) 

In the following we use the abbreviations L 2 = L 2 (D), L°° = L°°(D), 
Hq = Hq(D), H k = H k (D), < k < oo, for the standard Sobolev spaces. 
Let < ■,■ > and || • || = || ■ H2 denote the standard scalar product and 
norm in L 2 , respectively. Moreover, the norms for Hq, L°° are denoted by 
II " ||fr* 5 II " || 00 5 respectively. Due to the Poincare inequality (|T2|, p. 164), 
\\Atp\\ is an equivalent norm for Hq. It is well-known that the operator 
A = vA : L 2 -> L 2 with domain D(A) = H 2 (1 Hi is self-adjoint. Note 
that A generates an analytic semigroup S(t) on L 2 ([]19[). The spectrum 
of A consists of eigenvalues > Ai > A2 > A3 > . . . with corresponding 
normalized eigenf unctions </?i, <y9 2 , .... The set of these eigenfunctions is 
complete in L 2 . For example, for the square domain D = (0, 1) x (0, 1) the 
eigenvalues are given by —u(m 2 + n 2 )n 2 for m, n G N, and the associated 
eigenfunctions are suitable multiples of sin(m7rz) sm(mry). 

Now we can define an appropriate class of Wiener processes W. Let 
/3k{t), k G N, denote a family of independent real- valued Brownian motions. 
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Furthermore, choose positive constants p, k , k G N, such that 



rm — < 00 

k=l 1 1 

for some < 7 < 1. Then we consider the Wiener process W defined by 

00 

W(t) := <MW)<Pk , t > . (8) 
k=i 

We further assume that 

K (D) = inf inf ™as(DnB(x,y t p)) > Q> 

0<p<diam(D) (x,y)eD p 2 

where diam(D) is the diameter of .D (the least upper bound of two-point 
distances in D), meas(-) denotes the Lebesgue measure, and B(x, y; p) is the 
open disk centered at (x, y) and with radius p. We also assume that the 
eigenfunctions tp k satisfy 

if k E C (D), \<p k (x,y)\ < C, 

\d x tp k (x,y)\, \d y (p k (x,y)\ <Cy/\X^\, 

for (x, y) G D,k G N, and some constant C > 0. For the square domain 
D = (0,1) x (0,1), these conditions are all satisfied. Then, according to 
Theorem 5.2.9 in 0, the stochastic convolution 

W A (t) = [ S{t- s)dW(s) , t > , (9) 
Jo 

has a continuous version with values in Cq(D), the Banach space of contin- 
uous functions satisfying zero Dirichlet boundary conditions on D. 

If we define the nonlinear operator F by F{oS) = —ruj — (3ifj x — J(ifj,u), 
then (0) can be rewritten as the abstract evolution equation together with 
initial condition 

du = (Au + F(uj))dt + dW , (10) 
^(0)=a;o (11) 
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or in the mild (integral) form 

w(t) = S(t)w + I S{t- s)F{u{s))ds + W A {t) , (12) 
Jo 

where the stochastic convolution W A (t) is defined in (Q). 

By defining U(t) = uj(t) — W A (t), we obtain a deterministic mild (integral) 
equation 

U(t) = S{t)u + [ S(t- s)F{U{s) + W A (s))ds , (13) 
Jo 

or in its differential form 

U'(t) = AU(t)+F(U(t) + W A (t)) , (14) 
U(0)=u . (15) 

In the following we prove the local existence of U(t). We follow the 
approach in [[| or p. 261. We first show that the integral in ([HJ) makes 
sense for U G C([0, T]; L 2 ). Then, we obtain local existence for (|l"3"D by the 
Banach contraction mapping principle in L 2 . 

Note that since A generates an analytic semigroup S(t) on L 2 and has 
only negative eigenvalues, we have ([0, p. 74), for a > 0, 

S(t)(-A) a = (-AyS(t), (16) 

\\(-Ays(t)u\\ < £.|M|, (i7) 

||5(t)u|| < c-|H|. (18) 

Here and hereafter we use c to denote various constants. 

We first show that J * S(t — s)F(U(s) + Wa{s))cIs makes sense for U(-) + 
W A {-) (and thus [/(■)) in C([0,T];L 2 ). Recalling that u = U + W A , this 
follows from the following lemma. 

Lemma 1 Define the mapping T : C([0,T]; Hi) -> C([0,T]; L 2 ) &y 

{F{uj)){t):= I S(t- s)F{u{s))ds , tG[0,T], w € C([0, T]; L 2 ) . 

T/ien JF zs continuous, and it can be extended to a continuous mapping from 
the space C([0,T];L 2 ) to C([0,T];L 2 ). Furthermore, the image of the ex- 
tended mapping T is contained in C([0, T], H a (D)) for < a < |. 
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Proof: The continuity of T : C([0, T]; Hq) -> C([0,T];L 2 ) is obvious. As 
for extending the domain of T let u,u G C([0, T];L 2 ) be arbitrary. Using 
the abbreviations u> = Aip and u; = Atp we get 

F(a>) - F(u) = r(u - u) + (3(i> - tp) x + i> x (u - u) y + - ip) x u> y 

-4>y(u} - Uj) x + (i> - i>)yU x . 

Let a E [0,1), and consider an arbitrary (p £ D((—A) a ). Then the above 
identity implies 

/ = (i-A) a ip, S(t- s)F(u(s))ds - 1^ S(t- s)F(Q(s))ds^ 




(S(t - s)(-A) a <p, [r(u> + (3$ - ij,) x + ij x (u - u) y 



+ {i> - lfl)xUy - 1py(u - Uj) x + - 1p)yU x ](s))dS 

= [ (h + h + h + h + h + I^ds, 
Jo 

where 



h 


= <S(t- 


«)(- 






-u)](s) >, 


h 


= <S(t- 


«)(- 






- ^)x](s) >, 


h 


= <S{t- 


«)(- 


-A) a <p, 




- >, 


h 


= <S{t- 


*)(- 


-A)*<p, 


U- 


1p) X U)y\(s) >, 


h 


= <S{t- 


s)(- 






Q-u) x ](s) > 


h 


= <S(t- 


s)(- 






$)yU) x ](s) > . 



Now we estimate \Ik\, k = 1, . . . , 6, one by one, thereby omitting the argu- 
ment s. 

1^1 = \<S(t-s)(-A) a <p,r(Q-u)>\ 

< r\\S(t - s)(-A) a <p\\ ■ \\u - u\\ 

< rc(t - s)' a ■ \\<p\\ ■ 

|/ 2 | = \<S{t-s)(-A) a <p,l3(j>-1>) x >\ 

< 0\\S(t-8)(-A)*<p\\-\\$-il;) x \\ 

< (3c{t-s)- a - \\ip\\ • ||a>-a;||, 
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where we have used the Poincare inequality ([0, p. 164) on (jp — ip) x which 
has zero mean. Using the Cauchy-Schwarz inequality | < u, v > \ < \\u\\\\v\\ 
we also obtain 

|J 3 | = | < S(t-s)(-A) a ip,i> x (u-u) y > | 
= \<D v [(S(t-a)(-A) a <p)1> x ],Q-u>\ 

< | < D y (S(t- s)(-A) a ip)ip x ,u; - uj > | 

+ | < (S(t - s){-A) a p)ij X y,u - u > | 

< \\D y (S(t-a)(-A) a <p)iJ; x \\-\\u-u\\ 

+ \\S(t - s)(-A)VI|oo • H^ll • ||w - 
As for estimating 11^(5^ — s)(—A) a cp)ip x \\ we get 

\\D y {S{t-s){-A) a ^ x \\ < \\D y (S(t - s)(-Ay<p)h 4 ■ \\ 

< c \\D y (S(t-s)(-Ayip)\\ H i 

< c \\S(t-s)(-Ayip\\ s.\\oj\ 



4 

mm 



< c\\A^ +p S{t-s){-A) c 



< c{t-s)- ( 'i + 



■p+a) 



\UJ\ 



where we have used the inequality \\uv\\ < ||w||i||f||i, the continuity of the 

4 4 

3 1 reman 3 3 

mapping D y : — > ([Q, p. 56), the embedding D{{— A)* +p ) iJa for 
arbitrary p > p. 243), and the facts that Hz and if 1 are embedded 

in L 4 p. 217). Furthermore, 

||s(t- s )(-A)viloo < c\\s(t-s)(-Ayip\\ m+P 

< c\\{-A)i + >S(t-8){-AF<p\\ 

< c{t - s)-^ +p+a) ■ 



due to the smoothing property of the semigroup S and the embeddings 
D((-A)^ + f) ^ H 1+p ^ L°° for arbitrary p > (0, pp. 208, 243). Alto- 
gether we get 

|/ 3 | < [c(t " S y^ +p+a) + C(t - s)-(3+AH-«)] . U^ll . II^U . ||^ _ 

Similarly, the following estimates can be obtained: 

\h\ < [ c (t - s)-^^ + c(t - s)-^ +p+a) ] ■ H| • ||w|| ■ \\w-u\\ , 
|/ 5 | < [c(t - s)- ( f +p+a) + c(i - s)- ( ^+ a )] ■ II^H • \\u\\ ■ \\u-u>\\ , 
liel < [c(t - s)-^ +p+a) + c(t - s)-^ +p+a) ] ■ \\ip\\ ■ \\Q\\ ■ \\u-u\\ . 
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Thus, we have 

\I\ < [ (\h\ + \h\ + \h\ + \h\ + \h\ + \h\)ds 
Jo 

/ rc + pc , a 



1 - a 
+ 



8c 



sup \\co(s) — uj(s) 

0<s<t 

Ac 



l-Ap-Aa 



-p—a 



+ 



1 - 2p - 2a 



-p—a 



■\\ip\\- sup (\\uj(s)\\ + \\u(s)\\) ■ sup \\u(s) - u(s)\\, 

0<s<t 0<s<t 

provided the positive constants a and p satisfy < p + a < \. This finally 
implies 

/ S(t-s)F(u)(s))ds- [ S(t-s)F(u(s))dse D((-A) a ) for 0<a<], 

Jo Jo 4 



and 



< 



-A) a Qf S{t- s)F{u{s))ds - S(t - s)F(u(s))ds 
rc + f3c 1 



t a ■ sup \\uj(s) — u>(s 
1 — a o< s <t 



+ 



8c 



l-Ap-Aa 



Ac 



—p—a 



sup {\\u{s)\\ + \\u(s) 

0<s<t 



\-2p-2a 
sup \\co(s) — u>(s) || . 

0<s<t 



Especially for a = we obtain 



S(t-s)F(u(s))ds- / S(t-s)F(Q(s))ds 



< {rc + f3c)-t- sup ||u;(s) — ui(s) 

0<s<t 



+ 



Ac 



1-Ap l-2p 
sup (||a;(s)|| + ||a>(s)||) • sup ||a;(s) - a>(s) 

0<s<t 0<s<t 



for every < p < \. This completes the proof of the lemma. 



□ 
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We conclude from the above lemma that f* S(t — s)F(U(s) + Wa(s))cIs, 
considered as a mapping with argument U(-), can be extended to a bounded 
map from C([0, 77], L 2 (D)) into itself. 

Now we can follow p. 196 or 0, p. 201, to obtain that (|13[) has 
a unique local solution U(t), or (|10|) , (]TT|) has a unique local solution 
u(x,y,t), on [0,r), by the Banach contraction mapping principle. The solu- 
tion u(x,y,t) is in C([0, r]; L 2 (D)), as well as in C((0, r]; H a (D)), for arbi- 
trary < a < ~. 



3 Global solution processes 

In this section, we show that the solution U(t) is a priori bounded, in L 2 (D)- 
norm, on any finite interval [0,T]. This implies that the local solution U(t), 
and thus cu(x,y,t) is actually global in time. 

We consider (fil]) , fllEf) with replaced by a regular function V from 
the space C([0,T]; 77q(77)) and cj in 77(A) 

?/'(*) = Atf(t) + F(C/(t) + V(t)) , (19) 
U(0) = uq, 

where, we denote, U = Au, V = Av. More specifically, (0) is 

U' = vAU - r(U + V)~ (3(u + v) x - J(u + v,U + V). (20) 

Due to the smoothing effect of the sectorial operator A, and the fact that F 
is locally Lipschitz in U from H™ +1 n H 2m+2 to H™ n H 2m for m = 0, 1, 2, we 
conclude that the solution U of ( |T9"D is in i^Q fl i/ 2fc for fc = 0, 1, 2, and hence 
{7 is a strong solution (if V is smoother, the solution U is also smoother); 
see |IJJ, p. 73. 



We now estimate the norm ||J7(i)||. 

Multiplying ( PD by [7 and integrating over 7J, we get 



' (/ = -v [ \VU\ 2 -r [ (U + V)U 
Jd J d 



2dt 



-/3 / (u x + v x )U- / J(u + v,U + V)U 
Jd Jd 



■ [ |W| 2 -r f (U 2 + UV) 
Jd Jd 
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D 



-(3 / (u x U + v x U)- / J(u,V)U- / J(v,V)U 



D 



D 



v\\VU\\ 2 -r / (U 2 + UV)-[3 / (u x U + v x U) 



D 



D 



+ / (-U X VyU + UyV X U - V X VyU + VyV X U) , 

Jd 



(21) 



where we have used the fact that J D J(u, U)U = j D J(v, U)U = via inte- 
gration by parts; see also Hl7| . We estimate the right hand side of (|21~1) term 
by term. 



r / {U + V)U < r(l + c||y|| 00 )||[/|| : 



(22) 



D 



-(3 J (u x U + v x U) < P J \[ul + U 2 + v 2 x + U 2 } 

< Pc(\\u\\ 2 + \\v\\ 2 ) 

< (3c(\\U\\ 2 +\\V\\ 2 J, (23) 
where we have used the Poincare inequality on u x , v x , which have zero mean. 

/ (-U X VyU) = [ (U x U)yV 

Jd Jd 

= / (u xy UV + u x U y V) 
Jd 

< liviioo / \u xy u\+ [ (K| • ||V||oc)|c/,| 

Jd Jd 

< IIVIL / 1(U 2 +U 2 )+ [ (^\\V\\L + ^U? 



< ^ll^lloo^ + ill^lU^llf/f + ^llVt/ll 2 , (24) 

since § D u 2 xy is bounded by c f D (Au) 2 = c f D U 2 . We have also used the 
Young inequality (|[Z3|) to get that (\u x \ • WV^Uyl < i^lMIL + Wl for 
any positive real number e > 0. Similarly, we have 

D UyV x U < lll^Hoo (l + lllVHoo) W + |||W|| 2 , 
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-V X VyU) 



' D 



V X U)yV 

[ (v xy UV + v x U y V) 
Jd 



< 

< 
< 



VyV X U < 



D 



+ c 



D 



Y/ x \\v\\l + ^ y 



l\\v\\ 2 + 



+ c 

Putting (0) - (^Sj) into (fZT|), we get 
2 < {-v + 2e)-\\VU\\ 2 



3 

oo 
3 

oo 



4 

oo 



4 

oo 



(25) 



1 d 



+ [r(l + c 



+ /3c + c 



+Pc\\V\\l + 2c\\V\t + -c 



»(1 



l 

oo ■ 



+ 



Taking e = f , we finally obtain 



eft 



||f/(t)|| 2 < A(t)-||t/(t)|| 2 + J B(t), 



where 
A(t) 
B(t) 



r(l + c 



+ /3c + c 



2^11^11^+4011^111 



— c 



V 



> 0. 



+ 



(26) 

(27) 

>0, 
(28) 



Hence by the Gronwall inequality (f23j), we obtain 

t-t 

\\U(t)f < \\tu \\ 2 eti A{s)ds + / B{s)e£ A{T)dr ds, < t < T. (29) 

Jo 

Note that Hq(D) is embedded in Cq(D), the trajectories of WA_{t) can be 
uniformly approximated, on any finite interval [0, T], by functions V in 
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C([0, T}; Hq (D)), and D(A) is dense in L 2 (D). Thus the boundedness es- 
timate (p9|) is true for any local solution U(t) of (|14"D. This shows that the 
unique (local) solution does not blow up on any finite intervals. 
We thus have the following theorem. 

Theorem 1 For every initial condition uo(x,y) G L 2 (D), there exists a 
unique global mild solution u(x,y,t) of the quasigeostrophic model (fjj, 
and flfy. This solution is contained in the space C([0, T]; L 2 (D)) for 
every T > 0, as well as in C((0, T]; H a (D)) for allQ<a< \ and T > 0. 



4 Discussions 

There has been recent work on geophysical problems modeled by the ran- 
domly forced quasigeostrophic equation (e.g., PT| , [|T3]], JT5| , JT8] ], ]TD|). Un- 
like the deterministic quasigeostrophic equation whose well-posedness is well- 
known (e.g., [|J, p2| , Q), there seems no rigorous results on global existence 



and uniqueness of the randomly forced quasigeostrophic equation. Our work 
provides such a rigorous result on global existence and uniqueness, under 
very mild conditions. 
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